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ABSTRACT 

We investigate the possible origin of the strong radial gradient in the 
ultraviolet-to-infrared ratio in the nearby late-type spiral galaxy NGC300, and 
emphasize the importance of local variations in the geometry of the interstellar 
medium, concluding that they cannot be neglected with respect to metallicity 
effects. This analysis is based upon a combination of maps obtained both from 
space using GALEX (1350 to 2750 A) and Spitzer (3 to 160 /mi), and from the 
ground (UBVRI, Ha and H/5). We select ionizing stellar clusters associated with 
H II regions of widely varying morphologies, and derive their fundamental pa- 
rameters from population synthesis fitting of their spectral energy distributions, 
carefully measured to eliminate local background emission accurately. From these 
fits, we conclude that the stellar extinction law is highly variable in the line of 
sight of young clusters of similar ages. In the particular model geometry that 
we consider most appropriate to the sampled regions, we checked that our find- 
ings are not significantly altered by the correct treatment of radiative transfer 
effects. The variations are systematic in nature, in the sense that extinction laws 
of the Milky Way or LMC type, with a prominent 2175 A absorption feature 
and a modest far-ultraviolet slope, are associated with compact H II regions (the 
compacity being quantified in two different ways), while clusters surrounded by 
diffuse H II regions follow extinction laws of the 30 Doradus or SMC type, with 
a weaker 2175 A feature and a steeper far-ultraviolet rise. The Calzetti starburst 
attenuation law, although most often degenerate with the 30 Doradus extinction 
law, overpredicts ionizing photon fluxes by large amounts. We also find that the 
extinction law variations are correlated with the column density of dust species 
emitting in the near- and mid-infrared. Finally, we briefly discuss the nebular to 
stellar extinction ratios, and the excitation of aromatic band carriers, invalidating 
their claimed association with cold dust. 
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Subject headings: dust, extinction — galaxies: individual (NGC300) — galaxies: 
ISM — galaxies: star clusters — HII regions — ultraviolet: galaxies 

1. Introduction 

NGC300 is a gas-rich galaxy of Sd type, at a distance of only 2.1 Mpc (Freedman et al. 
1992). Its total extent in the blue light is over 20', and its inclination on the line of sight, 
derived from H I kinematics, is about 50° (Puche et al. 1990). It contains numerous bright 
H II complexes (Sersic 1966; Deharveng et al. 1988) and forms stars at a current total rate 
of about O.15M yr~ 1 (Helou et al. 2004). NGC300 was targetted for surveys of supernova 
remnants ; between 20 and 40 candidates were found by various selection methods (Blair & 
Long 1997; Pannuti et al. 2000; Payne et al. 2004), evidence that star formation has been 
going on for an extended period at a significant rate. Large numbers of Wolf-Rayet stars 
were also found (Breysacher et al. 1997; Schild et al. 2003). The oxygen abundance in the 
nebular gas phase varies between about 1.4 to 0.4 times the solar neighbourhood value, from 
the center up to a radial distance of about 10' (Deharveng et al. 1988). 

Many bubble-like structures are seen in the Ha emission, especially in the outer disk, 
where they reach diameters of ~ 600 pc. About 50% of the total Ha flux arises from the 
diffuse ionized gas component (Hoopes et al. 1996). Judging from the morphology of H II 
regions and the position of the most likely exciting sources, the hydrogen gas may be ionized 
very far away from young clusters (in projected distance). This argues for the existence 
of aged star-forming sites and a high degree of porosity of the interstellar medium around 
them. The large number of resolved stellar clusters and the variety of their associated H II 
regions, observed at a resolution of a few tens of parsecs, make it possible to study in great 
detail general patterns in the extinction and dust heating within a single galaxy. 

Deriving spatially-resolved extinction in galaxies of different types, and determining how 
the extinction law is affected by properties of the interstellar medium such as metallicity 
and geometric configuration, is a first step toward assessing the validity of common recipes 
used to correct for extinction the global ultraviolet emission of galaxies. This knowledge is 
fundamental to obtaining good prescriptions to estimate star formation rates and ages from 
current and future ultraviolet and infrared surveys. 

In this paper, we analyse jointly images of the stellar light from the far-ultraviolet to the 
far-red, of the hydrogen recombination emission in the Ha and H/3 lines, and of the infrared 
emission from 3 to 160 /mi. The data are described in Section 2. After deriving large- 
scale infrared-to-ultraviolet ratio and extinction maps (Section 3), we focus on individual 
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star- forming complexes selected in Ha and in the ultraviolet, and constrain the fundamental 
properties of ionizing stellar clusters, to assess how critically the extinction law and excitation 
of mid-infrared emitting dust at these spatial scales depend on the local geometry of the 
interstellar medium (Section 4). The results are discussed topic by topic in Section 5, and 
summarized in Section 6. 



2. Data 

2.1. Ultraviolet images 

NGC300 was observed with the GALEX space telescope in the FUV and NUV bands 
at effective wavelengths of 1516 and 2267A (Martin et al. 2005; Morrissey et al. 2005). 
GALEX has a circular field of view of 1.2° in diameter, and observations are performed in a 
spiral dither pattern. The data from the photon-counting, multichannel-plate detectors were 
processed with the analysis pipeline operated by the Caltech Science Operations Center, 
which converts the photon lists into flux-intensity images and applies flat-fielding and first- 
order geometric distortion corrections. The absolute calibration is derived from ground test 
data which show good agreement, to within 10%, with observations of white dwarf standard 
stars (Morrissey et al. 2005). See also the GALEX Atlas of galaxies (Gil de Paz et al. 2005), 
where NGC 300 is included, for more information. 

The angular resolution is of the order of 5" (FWHM), but varies with source brightness 
and position across the field. In particular, the point spread function in the FUV band is 
slightly elliptical (elongated along the south-east to north-west direction) in some parts of 
the map. A first distortion correction between the NUV and FUV maps was estimated using 
Legendre polynomials of order 4 and cross-correlating the positions of about 900 stars and 
clusters. Distortion residuals and point spread function variations were handled separately 
for each individual sub-field of size ~ 2', by fitting nearby point sources with gaussians in 
the FUV, NUV and U bands. 

2.2. Infrared images 

Broadband infrared images were acquired with the Spitzer space telescope at effective 
wavelengths of 3.6, 4.5, 5.7, 8.0, 24, 71 and 156 /zm, and were first discussed by Helou et 
al. (2004). The data processing and mosaicing is also described in this paper. The angular 
resolution (full width at half maximum) of the mosaics is 2.5 to 3" between 3.6 and 8/zm, 
5.7" at 24 urn, 16" at 70 /an and 38" at 160 fim. 
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The mid- infrared bands record gradually varying contributions, as a function of wave- 
length, from stars and dust, the global 3.6 fim emission being dominated by photospheres, 
the global 8 /iin emission deriving mainly from aromatic band carriers, and longer-wavelength 
bands covering essentially the thermal continuum from very small to big dust grains. 

The photospheric component was subtracted from the 8 /itn map as explained in Helou 
et al. (2004). A model for intermediate to old stellar populations was derived from the Star- 
burst99 population synthesis model (Leitherer et al. 1999), and was found to be remarkably 
invariant for solar and half-solar metallicities, and a broad range of star formation histories 
(excluding very young populations). Assuming that the 3.6 /zm emission is entirely of stellar 
origin, this map was then used to normalize the photospheric spectrum at each point, to 
which the relevant filter transmission curves were applied. This procedure may be inaccurate 
locally, because of the contribution from young stars, but the 8 /iin emission is heavily dom- 
inated by dust in star-forming regions, so that neglecting young stellar populations should 
not introduce any significant error. 



2.3. Optical images 

Mosaics in the UBVI bands of 20.5' x 20.5' were obtained from S.C. Kim (Kim et al. 
2004, private communication). An R-band mosaic of 25' x 25' and Ha+[N II] and H/3 maps 
in narrowband filters (in a field of view of 8.85' x 8.85') were acquired at the Las Campanas 
du Pont 2.5 m telescope with the direct CCD camera. The angular resolution in all the 
broadband optical images, after the mosaicing and the geometric transformations applied to 
align them, is of the order of 1.7" to 2" (FWHM). 

We also made use of an Ha mosaic covering the same field as the R-band mosaic to 
measure the Ha flux of the outermost H II regions ; the residual background was subtracted 
locally in each subfield separately. The map acquired by Hoopes et al. (1996), which was 
kindly provided to us by Charles Hoopes, supplied a means of checking the Ha+[N II] 
flux calibration, by comparing aperture photometry measurements on several bright H II 
complexes in both maps ; fluxes from the map of Hoopes et al. (1996) are systematically 
lower by about 10%. 

The H/5 photometric accuracy, however, is much more uncertain, because this line has 
a much smaller equivalent width than the Ha line, and larger errors are introduced by 
continuum subtraction and correction for stellar absorption features. Because of these un- 
certainties, if no flux correction is applied, the Ha/H/3 decrement values are unphysical in the 
case B recombination assumption, and in disagreement with published values for individual 
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H II regions. In order to correct them and derive an estimate of the nebular extinction from 
Ha and H[3 photometry, we scaled the decrement in such a way that extinction values derived 
by Webster & Smith (1983) from spectro-photometry for 5 regions could be reproduced. We 
used an intrinsic Ha/H/3 ratio of 2.86, appropriate for electronic densities and temperatures 
< 10 3 cm~ 3 and ~ 10 4 K, respectively, and replaced the extinction law assumed by Webster 
& Smith (1983) with the Cardelli et al. (1989) law. This procedure requires that the raw 
Ha/H/9 flux ratio be underestimated by 30% on average. Applying the above correction, the 
derived global extinction (from total fluxes in the H/3 field of view) is A(Ha) = 0.5 mag. 



3. Large-scale UV/IR ratio and extinction 

The ultraviolet disk reaches much farther out than the infrared disk. To quantify this 
statement, we convolved all maps to the resolution of the 160 /xm map, and estimated the 
total infrared emission (TIR) from a linear combination of the 24, 70 and 160 /xm powers (Dale 
& Helou 2002). To reach the 160 /xm full width at half-maximum, the maps at wavelengths 
up to 24 /xm were first convolved with a model of the point spread function at 70 /xm 1 ; 
these and the 70 /xm map were then convolved with gaussians of the appropriate widths. 
We computed surface brightness profiles in the FUV band and in the TIR, averaged in 
elliptical annuli to account for the inclination of NGC300 (Fig. 1). The FUV band was 
preferred over the NUV band for this comparison, because while the ratio of the extinction 
in these two bands is close to 1 (smaller than 1.3 for various extinction laws), the intrinsic 
FUV/NUV ratio of young stellar populations is > 1, the FUV band avoids the variable 
2175 A feature of the interstellar extinction curve (see Section 5.4), and the FUV band 
is minimally contaminated by foreground stars. On average, the FUV/TIR ratio increases 
exponentially with deprojected distance from the center, with irregularities mainly due to the 
FUV emission being less smooth than the infrared emission. The FUV/TIR ratio increases 
by at least a factor 5, maybe up to a factor 10, between the center and the outermost regions 
where it is still measurable. Also shown in Figure 1 is the opposite of the average slope of 
the metallicity gradient as derived by Deharveng et al. (1988). The slope of the FUV/TIR 
gradient is very similar, although somewhat smaller. At any given distance from the center, 
however, the dispersion in FUV/TIR ratios among individual star-forming regions is at least 
a factor 2, whereas the standard deviation of abundances is only of 0.11 dex, which seems 
to indicate that metallicity effects are not solely responsible for the FUV/TIR gradient. 
The UV-brightest regions tend to be deficient in infrared emission compared with their 
fainter counterparts. From studies of large samples of galaxies, a clear decrease of the global 



which can be found on the website: http://ssc.spitzer.caltech.edu/mips/psf.html . 
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ultraviolet to infrared flux ratio was found as a function of a measure of the star formation 
rate (Wang & Heckman 1996), or as a function of the blue luminosity (Buat & Burgarella 
1998) or the stellar mass (Pierini & Moller 2003). Interpreting the ultraviolet to infrared 
ratio variations within resolved galaxies could help constrain the underlying causes of this 
behavior. 

A stellar extinction map of NGC300 can be derived from a simple energy balance 
calculation, using all the maps between 1516 A and 3.5 /xm to reconstruct the stellar emission. 
The emerging fluxes in the FUV, NUV, UBVRI and IRAC1 (3.5 fim) bands were interpolated 
by realistic stellar spectral energy distributions, using the actual filter transmission curves. 
Assuming several different extinction laws, the power absorbed in this range is set equal to 
the total power emitted in the infrared. The derived A(FUV) ranges between 0.1 and 0.6 
(A(V) between and 0.23). Uncertainties of 20% in the flux calibration at both 70 and 
160 yum could increase the A(FUV) extinction values by 0.05 mag. 

The large H II complexes including shell-like structures correspond to regions of lower 
extinction, consistent with the idea that they are aged star-forming sites having partly 
dispersed the molecular and dust material around them. The correspondance with gas 
column density cannot be tested because no map in a molecular transition is available. 
The H I map of Puche et al. (1990), with a beam size of 17", allows us to see partial 
correspondance with the H I column density, though, especially in the outer disk (Fig. 5) ; 
the fact that there is no robust correlation with the extinction map is expected because H I 
does not trace adequately the total gas mass, which is usually dominated by the molecular 
phase in the inner parts of the disk where star formation is taking place actively. 

The radial profile of the ratio of the total stellar power (integrated between 1000 A and 
3.5 jum) to the total infrared power is notably different from the FUV/TIR profile: it shows 
a sharp decline from the center up to the pseudo-ring of most active star formation, reaching 
a minimum around a radial distance of 3 to 4', then rises again to reach its central value at 
a radial distance close to 10'. 

Several effects can contribute to the large-scale FUV/TIR gradient. 
1) The most obvious possibility is the variations in metallicity, which are expected to cause 
direct variations in the dust to gas ratio. However, various lines of evidence suggest that 
other parameters play a more fundamental role. Although the global dust to gas ratio of 
galaxies seems correlated with the far- infrared surface density (Andreani et al. 1995), which 
traces both the column density of dust and the star formation rate surface density, the dust 
to H I mass ratio of irregular galaxies is highly variable but not correlated with metallicity 
(Hunter et al. 1989). Sauvage et al. (1990) also showed that the far-infrared to blue flux 
ratio of dwarf galaxies is not correlated with their oxygen abundance. Within NGC300, 
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the nebular gas phase metallicity gradient, as determined by Deharveng et al. (1988), is 
slightly steeper than the FUV/TIR gradient. Interestingly, Bot et al. (2004) found that the 
dust to gas ratio of the Small Magellanic Cloud is much lower than expected from a linear 
scaling with metallicity, although the situation may differ between star-forming regions and 
the diffuse interstellar medium. 

2) The large-scale variations of the FUV/TIR ratio may also be affected significantly by 
the cumulative effects of local variations in the geometry of the sources. Several factors can 
contribute to systematically increase the distance between stellar clusters and the absorbing 
dust, notably: a decrease in the volume density of all the phases of the interstellar medium ; 
the maturity of the star forming regions, i.e. the outer disk could contain slightly more 
evolved young stellar populations which have had time to disperse gas and dust away. Aged 
star formation complexes indeed seem to efficiently increase the porosity of the surrounding 
interstellar medium, because they are associated with diffuse, shell-like H II regions, and are 
affected by lower extinction. 

3) The spectral energy distribution of the heating radiation could change in such a way that 
the total dust heating efficiency decreases. The fraction of the total stellar power emitted 
in the FUV band varies from less than 2% in the central region of the galaxy to about 5% 
in the outermost part of the disk, at and beyond a radial distance of 10' (the whole SED 
being corrected for extinction). The non- ionizing part of the stellar radiation thus becomes on 
average more energetic in the outer disk (the old populations have a smaller scale length than 
the younger populations emitting in the ultraviolet), and the dust heating balance should 
shift progressively from big grains emitting in the far-infrared to smaller grains emitting in 
the mid-infrared, which could explain part of the far- infrared flux deficit. We indeed observe 
a slight increase of the F24/TIR fraction in the outer disk (by about 30-40% compared with 
the annulus of intense star formation). However, this increase is too modest to account for 
the FUV/TIR gradient, which is mirrored closely by a similar FUV/F 24 gradient. 

Figure 2 shows the radial profile of the FUV/NUV energy ratio, which decreases on 
average beyond a radial distance of ~ 10'. Such a decrease of the FUV/NUV ratio in 
the outermost part of the disk cannot be explained by extinction or metallicity effects. It 
thus seems reasonable to assume that a gradient exists in the age of the youngest stellar 
populations which radiate in the ultraviolet. The range of FUV/NUV colors observed at 
large scales corresponds to average stellar population ages between about 16 Myr at radial 
distances near 10', and ~ 300 Myr in the outermost disk. 

In order to assess the role played by the geometry of the interstellar medium in the 
FUV/TIR variations, and to distinguish between the causes mentioned above, we now turn to 
individual H II complexes showing a wide range of morphologies, presumably corresponding 
to very different geometries. We take advantage of the high spatial resolution of the data 
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from the UV to the near-infrared, and characterize individual stellar clusters or small groups 
of clusters. If each one consists of coeval populations, then it is possible to lift the star- 
formation history degeneracy that affects the large-scale emission. 

4. Small-scale properties of H II complexes 

4.1. Selection of H II regions/young stellar clusters 

We selected individual stellar clusters or groups of clusters associated with prominent 
H II regions, that are visible in all the optical bands and sufficiently bright in the FUV and 
NUV bands to allow accurate photometry. Because of the selection method, the sample is 
biased against highly obscured regions, for which the analysis described in Section 4.3 would 
not be possible. Bright 24 /im sources which are not detected in the ultraviolet are discussed 
in Section 5.7. We also avoided objects with very complicated structure (making the iden- 
tification of the clusters ionizing a particular HII region ambiguous), with contamination 
by foreground stars or blending with nearby regions. The sample covers a large range of 
morphologies, mid-infrared colors, and UV to infrared ratios. They were labelled according 
to the H II region which they excite, in the catalog of Deharveng et al. (1988) (but the letter 
sublabels follow a different convention). Detailed notes, pertaining both to the morphology 
and to the modelling results, can be found in Appendix A. Images of the cluster fields in 
the U band, in Ha and at 8 jum are shown in Figure 6. 

Table 1 lists, for each group of clusters, the initial diameter of the stellar photometric 
aperture (see Sect. 4.2) and some observables: the Ha flux (corrected for an average [N II] 
contribution to the total Ha+[N II] flux of 10%), the nebular extinction derived from the 
Ha/H/3 decrement, and the Ha/R ratio measured within the stellar photometric area, R 
referring to line-subtracted broadband flux (i.e. pure stellar emission). Ha equivalent widths 
can be obtained by multiplying Ha/R by 1262 A (the equivalent width of the R filter). We 
also measured the F24/ F$ ratio at the peak of the 24 /zm emission (which can be offset from 
both the exciting clusters and the peak of the Ha emission), within a resolution element (2.5 
times the 24 /zm FWHM, except for confused regions for which the aperture was reduced to 
the FWHM in order to avoid neighboring sources), after convolving the 8/xm map to the 
24 fim angular resolution. 

Table 1 also indicates the association of the H II regions with supernova remnants in the 
catalogs of Pannuti et al. (2000), Payne et al. (2004) and Blair & Long (1997). We fail to 
see any direct connection between the presence of a supernova remnant and the morphology 
of the H II region, as some of these are classified here as compact (see in Section 5.4 the 
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Table 1. List of selected clusters, named according to the associated H II region in the 
catalog of Deharveng et al. (1988). A # symbol after the name of an object indicates 
association with a supernova remnant in the radio/X-ray catalogs of Pannuti et al. (2000) 
and Payne et al. (2004), and a * symbol indicates association with a supernova remnant in 

the optical catalog of Blair & Long (1997). 



name 


RA 


DEC (J2000) 


D 1 a 


F(Ha) 


4(Ha) 


Ra/R 






(h m s) (° ' 


) 


(arcsec) 


(10- 18 Wm- 2 ) 


(mag) 




(peak) 


D6 


00 


54 16.5 


- 37 34 56 


18.0 


139. 




0.44 


3.32 


D7-8 


00 


54 17.1 


— 37 35 19 


33.6 


652. 




0.15 


1.93 


D24 * 


00 


54 28.7 


— 37 41 36 


24.0 


549. 


0.38 


0.38 


1.30 


D39 # 


00 


54 38.5 


- 37 41 44 


12.0 


172. 


0.82 


0.87 


1.76 


D40 # 


00 


54 38.5 


- 37 42 41 


14.4 


173. 


0.15 


0.28 


0.94 


D45 * # 


00 


54 40.5 


- 37 40 55 


10.8 


363. 


0.81 


0.41 


0.93 


D53a * # 


00 


54 42.9 


- 37 43 12 


6.0 


156. 


0.79 


1.11 


3.42 


D53b * # 


00 


54 43.5 


- 37 43 11 


8.4 


365. 


0.69 


0.80 


2.64 


D53c * # 


00 


54 42.9 


- 37 43 00 


10.8 


380. 


0.39 


0.11 


~ 1 


D56 


00 


54 44.3 


- 37 40 27 


9.6 


109. 


0.85 


0.70 


1.00 




00 


54 45.4 


- 37 38 47 


12.0 




62 




1.21 


D76a # 


00 


54 50.4 


- 37 40 30 


8.4 


123. 


0.88 


0.69 


0.97 


D76b 


00 


54 51.0 


- 37 40 27 


7.2 


93. 


0.59 


0.43 


1.31 


D77 # 


00 


54 50.4 


- 37 38 24 


12.0 


407. 


0.38 


0.47 


2.07 


D79 # 


00 


54 51.2 


- 37 38 26 


16.8 


323. 


0.55 


0.46 


1.76 


D84 # 


00 


54 51.8 


- 37 39 39 


12.0 


316. 


0.66 


0.90 


2.22 


D88-90 


00 


54 53.4 


- 37 43 47 


33.6 


435. 


0.51 


0.32 


1.36 


D115 


00 


55 02.6 


- 37 38 28 


36.0 


327. 


0.37 


0.04 


1.21 


D118 # 


00 


55 04.1 


- 37 42 51 


30.0 


753. 


0.65 


0.33 


2.27 


D119 # 


00 


55 03.9 


- 37 43 20 


26.4 


985. 


0.49 


0.47 


2.22 


D122 


00 


55 04.7 


- 37 40 58 


14.4 


149. 


1.06 


0.63 


1.06 


D127 


00 


55 07.6 


- 37 41 51 


10.8 


104. 


0.70 


0.59 


0.94 


D129 * 


00 


55 08.9 


- 37 39 31 


15.6 


109. 


0.58 


0.51 


0.99 


D137a # 


00 


55 13.0 


- 37 41 41 


13.2 


919. 


0.41 


0.27 


3.43 


D137b * 


00 


55 12.4 


- 37 41 26 


28.8 


938. 


0.57 


0.20 


1.21 


D147 


00 


55 24.7 


- 37 39 33 


26.4 


148. 




0.10 


1.04 


D159 * # 


00 


55 33.9 


- 37 43 15 


24.0 


196. 




0.31 





a This is the diameter encompassing the pixels selected for the measurement of the cluster fluxes (see 
Section 4.2 for details on the photometric procedure). 
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discussion about a quantification of the compactness of the H II regions, and the significance 
of this parameter in interpreting the results). Among these catalogued remnants, weak X-ray 
sources are found close to D40, D45, and D53 ; the radio sources are all very weak, and some 
have flat radio spectral indices, consistent with thermal emission. Such is the case for the 
regions D61, D77, D84. The influence of these candidate young supernova remnants on the 
surrounding interstellar medium is thus unclear, and the shell-like structures visible in the 
Ha emission may correspond to older events. 

4.2. Photometry 

Special care was devoted to refine the photometric method so as to obtain accurate 
spectral energy distributions for the young clusters. In order to remove the local back- 
ground from old populations, we devised a method appropriate to the combination of UV 
and optical data. We had to take into account the facts that the underlying emission can 
vary significantly on the scale of the considered regions, that its brightness relative to that 
of the clusters is variable from one band to the other, increasing with wavelength and be- 
ing very high in the R and I bands, and that confusion makes any simple aperture always 
contaminated by old stellar populations and foreground stars. We define two diameters, 
constraining the source to be contained within the smaller (Di) and estimating the back- 
ground within the larger (D 2 ). Using the U band, where the background and contamination 
by old populations are minimal, we then determine refined photometric "apertures" in the 
following way. The average brightness b and standard deviation a of the background in D 2 
are computed iteratively in a narrowing brightness interval (±2(7; where a\ is the standard 
deviation at each step) until convergence is achieved. Then, the source flux is summed over 
the pixels above b + 3cr within Di and the background averaged over the pixels below b + 3a 
within D 2 (not excluding D\). Then, we use the exact same pixels to compute the source 
flux and background in all the other optical bands, which have been previously aligned to 
within 0.2arcsec, and control the results visually to ensure that no visible foreground star 
contaminates the aperture, and that D\ and D 2 were appropriately chosen. We find that 
this method produces significantly better results than by measuring the background in a 
concentric annulus or another aperture centered on a nearby empty patch. For clusters that 
appear single in the optical bands, we obtain good fits with a single-population model, which 
is not the case when we use the alternative methods above. 

Since the point spread function of the GALEX images varies across the field of view and 
is slightly elliptical in the FUV band, and because of residual distortion between the FUV 
and NUV images, we had to refine the photometric technique for the GALEX bands. In 
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Table 2. Broadband UV-optical photometric measurements, including the emission line 
contributions to the U and R bands. The given error bars are the formal measurements 
uncertainties only, to which we added 10% for the flux calibration uncertainty. 



name: 


FUV 

(/*Jy) 


NUV 


U 


B 


V 


R 


I 


D6 


892(±1%) 


819(±1%) 


835(±2%) 


603(±2%) 


540(±2%) 


575(±2%) 


391 (±4%) 


D7-8 


2419(±1%) 


2052(±1%) 


1442(±2%) 


1090(±2%) 


834(±2%) 


795(±2%) 


611(±4%) 


D24 


1377(±1%) 


1315(±1%) 


1275(±1%) 


898(±1%) 


730(±2%) 


809(±2%) 


482(±4%) 


D39 


287(±2%) 


257(±2%) 


392(±2%) 


205(±3%) 


166(±4%) 


288(±5%) 


110(±11%) 


D40 


510(±1%) 


489(±1%) 


439(±2%) 


309(±3%) 


234(±4%) 


272(±4%) 


133(±11%) 


D45 


526(±1%) 


613(±1%) 


862(±2%) 


574(±2%) 


472(±2%) 


645(±3%) 


326(±5%) 


D53a 


246(±1%) 


196(±2%) 


322(±2%) 


153(±2%) 


141(±3%) 


258(±3%) 


131(±6%) 


D53b 


733(±1%) 


632(±1%) 


913(±1%) 


519(±1%) 


452(±2%) 


712(±2%) 


376(±3%) 


D53c 


1387(±1%) 


1390(±1%) 


1182(±1%) 


971(±1%) 


849(±1%) 


933(±2%) 


846(±3%) 


D56 


225(±2%) 


200(±2%) 


276(±2%) 


165(±4%) 


134(±6%) 


213(±8%) 


80(±18%) 


D61 


398(±1%) 


356(±1%) 


518(±2%) 


301 (±2%) 


242(±3%) 


403(±4%) 


181 (±9%) 


D76a 


310(±1%) 


276(±2%) 


338(±4%) 


203(±5%) 


148(±10%) 


246(±12%) 


145(±14%) 


D76b 


360(±1%) 


301(±1%) 


274(±3%) 


189(±5%) 


147(±8%) 


214(±11%) 


132(±19%) 


D77 


1202(±1%) 


1082(±1%) 


906(±3%) 


606(±4%) 


503(±4%) 


700(±5%) 


419(±11%) 


D79 


929(±1%) 


817(±1%) 


908(±2%) 


593(±2%) 


482(±3%) 


717(±3%) 


398(±5%) 


D84 


411(±2%) 


354(±2%) 


621(±2%) 


381 (±2%) 


341(±4%) 


535(±5%) 


275(±9%) 


D88-90 


1183(±1%) 


1252(±1%) 


1273(±1%) 


880(±2%) 


700(±2%) 


877(±3%) 


491 (±5%) 


D115 


1713(±1%) 


1766(±1%) 


1516(±1%) 


1529(±1%) 


1332(±2%) 


1329(±2%) 


1326(±3%) 


D118 


2594(±1%) 


2587(±1%) 


2294(±1%) 


1688(±2%) 


1372(±3%) 


1645(±3%) 


914(±5%) 


D119 


2887(±1%) 


2771(±1%) 


3007(±1%) 


1887(±1%) 


1401(±2%) 


1882(±2%) 


894(±6%) 


D122 


252(±2%) 


239(±2%) 


369(±3%) 


232(±4%) 


177(±6%) 


292(±8%) 


153(±15%) 


D127 


163(±2%) 


139(±2%) 


212(±3%) 


137(±5%) 


119(±7%) 


188(±8%) 


138(±8%) 


D129 


369(±2%) 


309(±2%) 


354(±2%) 


210(±3%) 


161(±3%) 


245(±4%) 


105(±10%) 


D137a 


4124(±1%) 


3784(±1%) 


2966(±1%) 


2323(±1%) 


1853(±1%) 


1900(±1%) 


1182(±2%) 


D137b 


3397(±1%) 


3516(±1%) 


3625(±1%) 


2913(±1%) 


2446(±1%) 


2744(±1%) 


2046(±2%) 


D147 


925(±1%) 


835(±1%) 


615(±1%) 


560(±1%) 


433(±2%) 


410(±2%) 


276(±7%) 


D159 


696(±1%) 


664(±1%) 


554(±2%) 


354(±3%) 


296(±5%) 


406(±5%) 


228(±11%) 
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each small field around a selected cluster, we use point sources to fit positional offsets with 
respect to the optical images, and we determine the FWHM of the cluster independently in 
the FUV and NUV bands. In a few instances, where the selected cluster is partially confused 
by nearby clusters, we fitted and removed the latter before performing aperture photometry 
on the target. In a few cases when we have to truncate the aperture to avoid confusion, we 
apply an aperture correction, computed assuming that the source is perfectly gaussian. The 
measured broadband fluxes are listed in Table 2. 



4.3. Stellar population fits 

We used the Starburst99 synthesis model (Leitherer et al. 1999) to compute spectra 
of instantaneously-formed stellar clusters, with the metallicity Z = 0.008, well suited to 
NGC300, and a Salpeter initial mass function between 0.1 and 120 M . 

The transmission curves of the filters used in the observations were applied to the model 
spectra before doing any comparison between the models and the measured fluxes. For the 
optical filters, we used the definitions of the Johnson-Cousins system given by Fukugita et 
al. (1995) and the references therein. For the ultraviolet and infrared filters, we used the 
actual transmission curves of GALEX 2 and Spitzer 3 . 

We fit the observed FUV-NUV-UBVRI flux densities using three alternate extinction 
laws - Cardelli et al. (1989) for the Milky Way, with the ratio of total to selective extinction 
i?v = A(V)/E(B — V) = 3.1, appropriate for the diffuse interstellar medium ; Fitzpatrick 
(1986) for the LMC average and for the 30Doradus complex - and the attenuation law 
of Calzetti et al. (1994). The results should be interpreted while keeping in mind that 
the first three extinction curves were derived toward the line of sight of individual stars, 
and the last one is in principle only valid for the integrated emission of starburst galaxies, 
while we are applying them to spatial extractions within regions ranging between 60 and 
300 pc in diameter. While Misselt et al. (1999) have shown that the type of extinction law 
attributed by Fitzpatrick (1986) to 30 Doradus is in fact mostly observed in a region adjacent 
to 30 Doradus, the LMC2 giant shell, we retain the terminology of Fitzpatrick (1986) for 
reference to the particular extinction curve that he derived. Updated extinction curves for 
the LMC have been published by Misselt et al. (1999). Given that we are interested in 
extracting qualitative variations in the apparent local extinction law (see Section 5.4), and 



2 availablc from the website: http://galexgi.gsfc.nasa.gov/tools/Resolution-Response/index.html . 

3 available from: http://ssc. caltech. edu/irac/ 'spectraLresponse.html 
and http://ssc.caltech.edu/mips/spectraLresponse.html . 
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that the discrete extinction curves used here are in fact part of a continuum in dust properties 
(Gordon et al. 2003), these curves are sufficient for our purposes. 

We fit the cluster spectral energy distributions in a simple dust screen configuration, 
which is implicit in the four extinction curves considered here, i.e. ignoring the effects of 
scattering and radiative transfer through layers of dust mixed with the stars. These effects 
are extensively described by e.g. Witt et al. (1992) and Calzetti (2001). Allowing the relative 
geometry of stars and dust to be unconstrained would introduce insurmountable degeneracy 
with the fitted parameters of the stellar populations. But there is another motivation for 
assuming that the stars are not mixed with the dust. Our photometric technique, in effect, 
is designed to extract the emission from discrete stellar clusters and remove efficiently back- 
ground or foreground extended emission (Section 4.2). We argue that the stellar clusters 
can thus be treated as a collection of point sources, each with its own dust shell. This is 
especially true of the class of clusters that we call compact (see Section 5.4 for a quantita- 
tive definition). The smallest full width at half maximum of the images represents a linear 
distance of 20 pc, to be compared with the full widths at half maximum of young clusters 
in spiral galaxies as quantified by Larsen (2004), up to 10 pc. Even though the dust screen 
assumption does not reflect the likely geometry of the sources, we will show in Section 5.4 
that it is a good approximation for shell-type sources, for which the inclusion of radiative 
transfer effects (in particular scattering effects) does not alter our conclusions. Furthermore, 
the results concerning the variations of the attenuation law (which includes radiative transfer 
effects, as opposed to the pure extinction law) are independent of any assumptions about 
the source geometry. Note that scattering should be a negligible component of the Calzetti 
et al. (1994) law, because light scattered in the line of sight compensates on average light 
scattered out of the line of sight for regions of large size (several kiloparsec). 

The flux calibration uncertainty is assumed to be 10% for all the filters, to which the 
formal measurement uncertainty is added. The fits were not constrained to reproduce the 3- 
5 /iin data, but only not to overproduce them (within the error bars) , even when these bands 
are dominated by stellar emission. We find in some cases that no single coeval population 
can reproduce the photometric data, in which case we allow two populations of different ages 
but same extinction. This occurs for the double cluster D53c and for some complex fields 
associated with very large H II regions (D137b, D115). All the other regions were fitted with 
a single population. 

In H II regions, the Ha and [O II] lines at 6563 A and 3727 A, respectively, contribute 
a large fraction of the total flux in the R and U bands (up to 50% in R, and up to 40% 
in U, respectively, within the chosen apertures), except in very diffuse H II regions. The 
contribution of the H/3 and [O III] lines can be neglected because they are observed through 
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low-transmission parts of the B and V filters, and are intrinsically weaker. We accounted for 
Ha and [O II] fluxes by subtracting scaled versions of the Ha map from the R and U maps, 
using the [O II] /Ha ratios measured by Webster & Smith (1983) and d'Odorico et al. (1983) 
and correcting them for differential extinction. Whenever a measurement of the [O II] /Ha 
ratio is lacking, we adopt the average ratio of the H II regions discussed by Webster & Smith 
(1983) and d'Odorico et al. (1983) ([O II]/Ha = 1.13). Because the [O II] contribution 
estimated in this way is uncertain and because of possible aperture mismatch, we affect a 
smaller weight to the U-band data than to the other bands, recovering an updated [O II] /Ha 
ratio from the SED fits (Table 3). 

The output parameters of the fits are the age(s) of the clusters, their mass(es), and 
the monochromatic extinction at 1516 A, the effective wavelength of the FUV filter. The 
extinction and ages take discrete values, varying in steps of 0.1 mag and lMyr, respectively. 
For each extinction, each age of the young clusters, and optionally each age of older clusters, 
the masses have to be solved from a set of 7 equations (for the 7 bands from the FUV to 
I). In order not to bestow a larger weight on any particular combination of bands, we solve 
separately each equation with the new equation resulting from the sum of all the fluxes, 
and then we compute the average of the 7 solutions, with a lower weight for the U band. 
The ionizing photon flux predicted by the best fit is compared with the ionizing photon flux 
estimated from the Ha flux, corrected for extinction using the Ha/H/3 decrement, but is left 
unconstrained by the observations. We usually find agreement to within 20%, with a few 
exceptions (Table 3). To determine the confidence intervals of the cluster parameters, the 
various fits are only constrained to reproduce the ionizing photon flux from the best fit to 
within 30% (we set a constraint on the error, but not on the absolute value). 

We expect the error on the observed ionizing photon flux to be large, because of un- 
certainties on the Ha/H/3 decrement, because of uncertainties in the extent of the region 
directly under the influence of the ionizing stars, and because an unknown fraction of the 
ionizing photons may leak out of the H II regions. The few cases where the agreement be- 
tween the model value and the observed value is poorer than 20% may be due to local errors 
in the subtraction of the continuum and stellar absorption from the B.(3 line, or to increased 
porosity of the surrounding interstellar medium. One of the regions for which the best-fit 
model overpredicts the measured ionizing photon flux by a large percentage is D137b (by 
40% for the 30Doradus law, and 100% for the Calzetti et al. (1994) law). This region also 
stands out as the most peculiar in terms of morphology of the hydrogen gas, suggesting that 
previous generations of stars in the recent past have injected lots of mechanical energy in 
the surrounding medium, carving a tunnel-shaped void. It can thus be hypothesized that a 
sizeable fraction of the ionizing photons are leaking out of the visible H II region. According 
to studies of the diffuse ionized gas in spiral galaxies, up to 50% of the ionizing photons may 
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Table 3. Results of the population synthesis fitting. 
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a When no error bar on the age is given, it means that the la interval is smaller than the time step of 
1 Myr used in the fits. 

b Since the age and mass of secondary stellar populations are always ill-constrained, this table indicates 
only the best-fit values. 

c The first value is the input (from Webster & Smith (1983) or d'Odorico et al. (1983), or the average ratio 
of the regions discussed in these papers when followed by a star) ; the second value is the output of the best 
fit (not given if the Ha/R ratio is below 0.15, too low to derive [O II] /Ha from the fits). All the sub-regions 
of D53 were assumed to have the same [O II] /Ha ratio. 

d The ratio of the ionizing photon flux predicted by the model to that estimated from the observations is 
given for the best fit of each applicable extinction law, in the following order: Milky Way, LMC, 30Doradus, 
and Calzetti et al. (1994). Whenever the 30Doradus and Calzetti et al. (1994) laws both provide acceptable 
fits, since the Calzetti et al. (1994) law predicts an ionizing photon flux usually in excess, the cluster 
parameters and their confidence intervals are derived only for the 30Doradus law (the Calzetti ct al. (1994) 
law produces larger masses and higher extinctions). 
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escape H II regions (Ferguson et al. 1996). Alternatively, the inclusion of a second population 
increases the degeneracy of the fits, and the fitted stellar extinction could be too high. 

We assess the validity of the different extinction laws by first comparing the x 2 values 
of the best fits. But since the absolute uncertainty was fixed arbitrarily to 10%, not taking 
into account the smaller relative uncertainty between the two GALEX bands, we consider 
that a x 2 value close to unity does not necessarily indicate that the fit is acceptable. We 
also check its ability to reproduce colors, in particular the FUV-NUV color, which is the 
most useful discriminant between the different extinction laws. The ratio of the extinction 
in these two bands varies from ~ 1.00 for the Milky Way and ~ 1.10 for the LMC to ~ 1.27 
for 30Doradus and ~ 1.20 for the Calzetti et al. (1994) law. 

Whenever some significant interstellar component at 3-5 /zm was indicated by the pop- 
ulation synthesis model, we confirmed this finding by examining the morphology of the dust 
and stellar clusters (i.e. by checking that the structure of the 3-5 /xm emission is much more 
similar to that at 8/im than to that of the stars). 

5. Results and discussion 
5.1. Ages 

One characteristic of the fits assuming instantaneous star formation is that they nearly 
always predict a cluster age of 3 Myr, with the notable exception of D53a, which is fitted 
with a stellar population of 1 to 2 Myr, and is one of the regions with the highest F 2 4/F 8 flux 
ratios. This uniformity of age is likely an artifact of adopting too simplified a representation 
of the clusters (instantaneous star formation model with coarse time steps of 1 Myr), whereas 
their formation may go on for several Myr instead of being coeval. That the clusters are all 
young and observed before the bulk of supernova explosions can be inferred from the fact 
that they produce a significant amount of ionizing radiation, since most were selected to 
be associated with bright H II regions. The case of D53a is interesting, as it confirms that 
complexes containing younger stars tend to have higher ratios of very small grain emission 
to aromatic band emission (Sauvage et al. 1990; Roussel et al. 2001b; Dale et al. 2004). 

5.2. Masses 

When considering such small stellar populations as those extracted here, formed of one 
or a few young clusters, statistical fluctuations of the mass spectrum of massive stars may 
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become severe enough to affect the results of the fits. For the adopted initial mass function, 
the ratio N(M > Mh ot )/M tot , where the numerator is the number of stars more massive than 
the threshold Mh t and the denominator is the total cluster mass as given in Table 3, is equal 
to 0.1264M^ 1 t 35 - 0.0002, with masses in units of solar mass. The whole range of cluster 
masses in Table 3 thus corresponds to between 9 and 126 O stars with Mhot = 16 M Q (Vacca 
1994). Adopting as a working definition of ionizing stars Mhot = 10 M (Kennicutt 1984), 
the clusters masses also correspond to between 17 and 245 ionizing stars. The sampling of 
the initial mass function is thus a concern for the smallest clusters analysed in this paper, in 
which a single ionizing star may account for up to 5% of the far-ultraviolet emission. Note 
that a few regions, for which we could not obtain any satisfactory fits, were discarded from 
our sample ; one of the possible reasons for this is small-number statistics of ionizing stars 
in low-mass stellar clusters. 



5.3. Stellar and nebular extinctions 

The fitted stellar extinction at the effective wavelength of the FUV filter varies between 
nearly zero to ~ 1.5 mag. Its correspondence with the nebular extinction is shown in Figure 9. 
The lack of a tight correlation between the two can be understood in terms of spatial offsets 
between H II regions and clusters, and highly variable geometric configuration from one 
region to another. Within this sample, the A(H«)/Ar ratio ranges from 1 to more than 3.5 . 
The deviation from the value usually adopted for actively star-forming galaxies (Calzetti et 
al. 1994) can thus be greater than ±50% for individual regions. 

There is a good qualitative correspondance between the extinction map at a resolution 
of 38" and the extinctions of individual regions derived from the population synthesis fits, 
considering that significant deviations are expected from the fact that extinction is highly 
variable on small scales. It should also be noted that while the population synthesis fits give 
the extinction of the young stellar clusters only, the extinction in Figure 4, by construction, 
includes all the stellar populations at any given location. Since intermediate-age and old 
populations are, on average, much less affected by extinction at any wavelength, because 
decoupled over time from their parent gas clouds, the extinction values derived from the 
global energy budget are expected to be much lower than those of the ionizing clusters ; this 
is indeed observed. 
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5.4. Extinction law and geometry of the interstellar medium 

We find clear differences in which extinction laws are applicable, as a function of the 
geometry of the H II regions with respect to the exciting clusters. The most compact regions 
follow either the Milky Way or the LMC law, whereas the clusters associated with H II 
regions forming shells and diffuse structures follow either the 30 Doradus or the Calzetti et 
al. (1994) law. Usually, the degeneracy between the latter two cannot be lifted, but models 
using the Calzetti et al. (1994) attenuation law systematically overpredict the ionizing photon 
flux by much larger amounts than the 30 Doradus law. We therefore consider the 30 Doradus 
law to yield a more reliable representation of the data, even though the stellar SED fits are 
formally of comparable quality using either of these two laws. 

To qualitatively ascertain the consequences of neglecting radiative transfer effects, we 
examine the predictions of the DIRTY model (Witt & Gordon 2000; Gordon et al. 2001; 
Misselt et al. 2001) with empirically-derived dust properties, in the configuration of a dust 
shell surrounding a stellar cluster, appropriate for our sample (Section 4.3). The dust prop- 
erties of the LMC2 giant shell (called 30 Doradus here), not included in the original models, 
were derived by Gordon et al. (2003) 4 . The results obtained from the pure extinction laws 
for the Milky Way, the LMC2 region and the SMC are compared with those obtained from 
the corresponding attenuation laws for two types of local dust distribution: homogeneous 
and clumpy, the latter using the parameters defined by Witt & Gordon (2000). While other 
geometric configurations are available within DIRTY, better applicable to the extended emis- 
sion of large regions within galaxies, it is beyond the scope of this observational study to 
perform comprehensive modelling. We stress that the basic result of large variations in the 
attenuation law, whether they correspond to genuine variations in the extinction law or not, 
is independent of any assumptions about the relative geometry of stars and dust. 

Figure 10 shows the extinction and attenuation curves in the model shell configuration, 
for a particular optical depth, which corresponds to A 1516 = 0.7 using the Cardelli extinction 
law, within the range of derived extinction values for our clusters. The inclusion of scattering 
produces emerging spectra that are bluer, with significantly higher FUV/NUV ratios, as more 
light is scattered in the line of sight than out. It is thus conceivable that the clusters which 
were well fitted with the Cardelli extinction law could be equally well fitted with a LMC or 
SMC-based attenuation law. However, even at the low optical depth considered in Figure 10, 
the attenuation curves applicable to the clumpy case are extremely grey, and thus would not 
provide any good fit to the observed spectral energy distributions, since they would predict 
too much intrinsic emission in the optical bands. 



The data tables were retrieved from the DIRTY website: http://dirty.as.arizona.edu/ . 
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We now examine the aptitude of each attenuation law to reproduce the observed FUV/NUV 
flux density ratios. Figure 11 shows the predictions from the DIRTY model. The intrinsic 
FUV/NUV ratio of a 3Myr-old cluster is 1.18, and it varies between 1.13 and 1.31 for ages 
between 1 and 5Myr. This uncertainty is reflected in the error bars in Figure 11. 
Cluters well fitted with the Cardelli or LMC extinction law: 

The data points overlap only with the Milky Way extinction law and the LMC2 extinction 
and attenuation laws. However, the attenuation law using the LMC2 dust properties, in 
the homogeneous case, is not much different from the pure extinction law, i.e. the effects 
of scattering are almost negligible. It is in principle possible to reproduce the observed 
FUV/NUV ratios with the LMC2 or 30Doradus laws, but only if the clusters are younger 
(1 or 2Myr). This seems unlikely, since the predicted ionizing photon fluxes would then be 
in strong excess with respect to the observational estimates (which are more robust for the 
compact HII regions associated with these clusters than for the diffuse HII regions). We note 
that the effects of scattering are strongest for the Milky Way dust properties, and that none 
of the clusters would be well represented in the dust shell configuration with the Milky Way 
attenuation laws (i.e. including scattering), either in the homogeneous or in the clumpy case. 
Cluters well fitted with the 30 Doradus extinction law: 

The data points lie in between the LMC2 (30 Doradus) and the SMC laws. Again, the SMC- 
dust attenuation law, in the homogeneous case, is very close to the pure extinction law, and 
the effects of scattering are also very small in this configuration. But the extinction laws in 
the line of sight of these clusters may span the full range between the 30 Doradus and the 
SMC laws. 

We conclude that the correct treatment of radiative transfer effects would not drastically 
alter our qualitative conclusions in what follows, unless stars had extended spatial distri- 
butions and were significantly mixed with the dust, which we argue is very unlikely for the 
stellar clusters that we classify as compact. 

To quantify the compactness of the regions in a convenient way, we use the Ha/R flux 
ratio, measured in the area defining the stellar aperture (i.e. if the H II region extends 
far away from the clusters, only a fraction of the Ha emission is considered in this ratio). 
While this is proportional to the Ha equivalent width, we argue below that it is also a 
valid measure of the compactness of the H II region. Although it is potentially distorted by 
differential extinction between the nebular phase and the stars, this effect should not have 
a large impact. Using the results from Figure 9, we derive that for the highest extinctions, 
observed values of Ha/R close to 1 at one extreme or close to 0.2 at the other extreme may 
be shifted by up to 0.9 upward or 0.2 upward, respectively, when corrected for extinction. 
Since the regions with the highest extinctions are also those with the highest Ha/R values, 
correcting rigorously for extinction would preserve the trend discussed below. Figure 12 
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confirms that the Ha/R ratio provides a reliable quantification of the degree of compactness 
of the entire H II region. We prefer to adopt the Ha/R indicator, because it is a relatively 
unbiased measurement: unlike the ordinate of Figure 12, it does not imply any assumption 
about the total extent of the H II region, which is difficult to define owing to the presence 
of neighboring H II regions and the diffuse and filamentary aspect of the Ha emission. 

Figure 13 shows the dependence of the extinction law on the Ha/R ratio. Whenever a 
region sits in between two extinction laws, it means that its SED is formally as well repre- 
sented with either law as with the other. Unconstrained regions are those whose SED is well 
fitted by at least three of the four extinction laws (usually because of very small extinction). 
To indicate cases when the measured ionizing flux is overestimated by a large amount by the 
Calzetti et al. (1994) attenuation law, the points are plotted closer to the 30Doradus law 
locus. The main discriminant between the various extinction laws is the A(FUV)/A(NUV) 
ratio (Sect. 4.3). The increase in A(FUV)/A(NUV), going from the Cardelli et al. (1989) 
law to the LMC law to the 30 Doradus law, is due to both the gradual disappearance of the 
so-called 2175 A bump, which is covered by the NUV band, and the steepening of the FUV 
rise. Although the Calzetti et al. (1994) law is devoid of the 2175 A bump, the correspond- 
ing A(FUV)/A(NUV) ratio is slightly smaller than that of the 30 Doradus law, because the 
Calzetti et al. (1994) law has a flatter wavelength dependence. 

That the extinction law should vary within a galaxy is not a complete surprise given 
that some significant variations have been observed within the Milky Way and the Magellanic 
Clouds. In the SMC, Lequeux et al. (1982) discussed the line of sight to a star with a Milky 
Way- type extinction law. In the LMC, Misselt et al. (1999) discussed a line of sight associated 
with a giant shell, with a much weaker 2175 A bump than its immediate surroundings. In 
the Milky Way, LMC-type extinction laws have be found toward low-density lines of sight 
(Clayton et al. 2000), and a SMC-type extinction law was derived toward a region suggested 
to have been profoundly affected by supernova shocks (Valencic et al. 2003). Although the 
disappearance of the 2175 A bump, together with the steepening of the far-ultraviolet rise, 
are usually ascribed to a decrease in metallicity, which could explain the differences between 
the average extinction laws of the Milky Way and Magellanic Clouds, Mas-Hesse & Kunth 
(1999) found no such dependence of the extinction curve on metallicity, for a sample of dwarf 
star-forming galaxies. As pointed out by Misselt et al. (1999), an apparent dependence on 
metallicity may be only an indirect effect, because the dust grain properties and geometry 
of the interstellar medium are expected to be on average related to the metallicity. 

It has been suggested that the intensity of star formation in the local environment 
could be the main driver of the extinction law shape, causing the destruction or processing 
of the dust grains responsible for the 2175 A feature, via radiation hardness and supernova 
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shocks (Fitzpatrick & Massa 1986; Gordon et al. 1997; Gordon & Clayton 1998; Mas-Hesse 
& Kunth 1999). This scenario is well in line with the attribution of part of the 2175 A feature 
to aromatic hydrocarbons (Joblin et al. 1992; Duley & Seahra 1998), as this dust species is 
known to be less resilient than very small grains and to be destroyed by intense radiation 
fields. Supporting this identification, Vermeij et al. (2002) have reported variations of the 
mid-infrared band ratios (caused by changes in molecular structure, or degree of ionization 
and hydrogenation) that are correlated with the strength of the 2175 A bump. It has been 
proposed that aromatic hydrocarbons are also responsible for part of the far-ultraviolet rise 
(Puget & Leger 1989), and very small carbonaceous grains (of three-dimensional structure) 
are sometimes assumed to be the main contributor to the 2175 A bump (Draine & Lee 1984; 
Hecht 1986; Sakata et al. 1994). 

The present work, however, supports the idea that geometry effects are much more im- 
portant than the distinction between quiescent and actively star-forming regions. According 
to the SED fits, all the exciting clusters have very similar ages. Therefore, rather than being 
a pure indicator of the evolution of the clusters, as the equivalent widths of the H recom- 
bination lines are usually assumed to be, the Hct/R ratio can be interpreted here as mainly 
governed by the geometric evolution of the whole region, bearing traces of previous genera- 
tions of stars in the shaping of the interstellar gas. Indeed, the giant shells of hydrogen that 
are observed in many locations across the disk of NGC300, where they are illuminated by 
ionizing clusters, may have been carved by supernova explosions, in which case they signal 
past star formation at the same location, more than 5Myr ago. The difference in geome- 
try between the regions D53a-b and D53c, which are part of the same complex and whose 
exciting clusters have approximately the same age, is particularly significant. High F 2 4/F 8 
ratios, indicative of intense radiation fields, can also be found both in compact and diffuse 
star formation complexes. The most obvious difference between clusters associated with a 
Milky Way-LMC extinction law and those associated with a 30 Doradus-Calzetti law is in 
the distribution of the material around them. For the latter, both the atomic hydrogen and 
the mid-infrared-emitting dust (aromatic hydrocarbons and very small grains) have been 
removed far away, presumably by stellar winds and supernova shocks from previous gener- 
ations of stars at the same location. This does not preclude important differences as well 
in the processing of dust grains by supernova shocks, but if the modified dust is dissipated 
together with the gas, it is only weakly participating in the extinction of the clusters. 

We can speculate that the effects of a varying extinction law within and among galaxies 
may resemble the effects of a varying initial mass function (Rosa & Benvenuti 1994), since 
for a given optical depth, the emerging ultraviolet continuum is much flatter in the case of 
a 30 Doradus-Calzetti law than for a Milky Way law. An upper mass cutoff of the initial 
mass function is not supported in NGC300, since Wolf-Rayet features have been observed 
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throughout its disk (Schild et al. 2003). 

Finally, we can expect the behavior of more massive galaxies and of starbursts, where 
H II regions may generally be more centrally distributed and more compact than in NGC 300, 
to be significantly different. 

5.5. Interstellar emission in the near-infrared 

In some regions, we detect a clear excess at 3-5 /im above the stellar model. The improve- 
ment with respect to the previous treatment in Helou et al. (2004) (where it was assumed 
that the 3.5 /zm emission covered by the IRAC1 band is entirely stellar) allows us to derive 
more accurately the amount of interstellar emission at these wavelengths. However, our 
method restricts the detectable interstellar component at 3-5 /im to the exact location of 
the stellar emission ; if the hot dust in this wavelength range is significantly offset from the 
exciting source, it will not be detected. The 3.5/zm/4.5/zm color of young clusters is very 
similar to that of dust (close to 1), making it very difficult to separate the two components if 
extensive information at other wavelengths does not exist. Complexes where the interstellar 
component at 3-5 fim is high with respect to the stellar fluxes (within the stellar aperture) 
are among the most compact (Fig. 14 and 13). 

To explore the link between the extinction law and the dust emission, we adopt mea- 
surements of the 3-5 fim interstellar emission, as the favorable angular resolution in these 
bands allows us to select the dust emission within the stellar "apertures" , directly in the line 
of sight of the clusters. Figure 14 confirms that switching from the Ha equivalent width to 
the column density of near-infrared-emitting dust in the line of sight of the clusters gives the 
same trend as seen in Figure 13: the shape of the extinction law seems to depend importantly 
on the amount of interstellar material which re-emits the absorbed energy in the near- to 
mid-infrared. We cannot explore related variations in the mid-infrared to far-infrared colors, 
because the angular resolution at far-infrared wavelengths is too coarse. 

Even though this 3-5 /im dust component may not be physically related to the other mid- 
infrared dust species, the aromatic band carriers and very small grains, Figure 15 shows that 
it scales reasonably well with the 8 /im aromatic band emission. These measurements suffer 
from the degraded angular resolution at 8 /im combined with the narrow stellar "apertures" 
used to fit and extract the interstellar component from the 3.6 /im band, but are nevertheless 
in line with the F 3 . 6 /F 8 flux ratios derived by Lu (2004) for interstellar emission in star- 
forming galaxies, as well as for the general interstellar medium of the Milky Way. 
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5.6. Porosity of the interstellar medium and excitation of aromatic band 

carriers 

Examination of the morphology of the H II complexes in Ha provides evidence that the 
interstellar medium is highly porous in some regions, thus that the 8 /xm-emitting dust can 
be heated by photons travelling as far as several hundred parsecs (since the NUV photons 
travel farther away than the ionizing photons, which are able to excite Ha emission a few 
hundred parsec away from the clusters). Diffuse emission also exists at 24 /im, but is harder 
to detect, because it is intrinsically fainter than at 8 /im and observed at much lower angular 
resolution, increasing the confusion by nearby bright sources, and the 8 /im observations are 
deeper. Although the 8 /im/Ha flux ratio varies dramatically from bright H II regions to 
more diffuse regions (Helou et al. 2004), and the 8 /im emission has a more filamentary aspect 
than the Ha and 24 /im emission, virtually every 8 /im peak is associated with a close-by Ha 
peak (see Figure 6). 

The claim of Haas et al. (2002) that aromatic band emission is more closely correlated 
with very cold dust than with hot dust within galaxies does not withstand examination 
in NGC 300. Figure 16 shows the relation between the aromatic band emission at 8 /im 
and the emission in the far-infrared, at 70 and 160 /im. As a comparison, the relation 
between 24 /im and the far-infrared emission is also included. The 160 /im emission can be 
considered a tracer of the cold dust, and the 70 /im emission arises from dust that is on 
average hotter, with a significant contribution from impulsively-heated very small grains 
(Desert et al. 1990). Figure 16 indicates clearly that aromatic bands are not associated 
with cold dust, and are not excited by the same radiation as the big grains. In the 8/im- 
70 /im diagram, at high surface brightnesses, a branch with a slightly lower slope than in 
the low-brightness part is visible. This branch is slightly more prominent in the 8 /im-24 /im 
diagram ; it is linked with the disappearance of the aromatic bands from inside the H II 
regions. Nevertheless, the 8 /xm-70 /im and 8 /iin-24 /iin relations are much closer to linearity 
than the 8 /iin-160 /im relation. The linearity between 8 /im and 850 /im surface brightnesses 
reported by Haas et al. (2002), obtained at low angular resolution in galaxies of moderate 
angular size, and more importantly at mediocre sensitivity at 850 /im, may simply be a 
consequence of extracting only the brightest and most active regions, which generally have 
dust temperature distributions shifted toward higher values, and discarding more diffuse 
regions, with cooler dust temperatures. Their result thus does not imply that the aromatic 
bands and the cold dust component are physically related, nor heated by the same radiation. 



-24- 



5.7. Embedded and bare sources 

Within the inner disk, a few non-stellar sources bright at 24 /zm have weak optical and 
no ultraviolet counterparts (Table 4). Two of these are probably background distant galax- 
ies, judging from their morphology in optical high-resolution images. The five other sources 
coincide with very faint Ha emission, have no other optical counterpart, and may be ultra- 
compact H II regions. We consider it likely at least for the object catalogued by Deharveng 
et al. (1988) as the faint and small H II region D125. If the Ha/H/3 decrement of D125 
can be trusted, it indicates a lower limit of the extinction in the Ha line of 2.8 mag, much 
higher than for the H II regions associated with the clusters examined here (see Table 1). 
Its mid-infrared counterpart has the highest F 2 4/F 8 ratio within the galaxy (~ 4.5) and is 
still visible at 70 /iin despite the degraded angular resolution. 

Conversely, some far-ultraviolet and Ha sources are not associated with any detectable 
mid-infrared dust emission. These regions are of low brightness and diffuse in the Ha emis- 
sion, suggesting that they are aged H II regions, excited by slightly evolved clusters (see the 
previous paragraph). Both types of sources contribute negligibly to the emission from the 
optically-bright disk, but the latter (bare sources) become predominant in the outer disk. 

6. Summary and conclusion 

The large-scale UV/IR ratio gradient observed in NGC300 is probably the consequence 
of a combination of factors, all causing an increase in the mean free path of ultraviolet 
photons. In this paper, we have focussed on geometric factors, and shown that they cannot 
be neglected. 

Throughout the disk of NGC300, mid-infrared peaks are commonly offset from both 
the ionized gas peaks and stellar clusters, except for compact regions (Helou et al. 2004). 
This is because we are able to partially resolve H II regions from photodissociation regions. 
However, in the outer disk, even the far-infrared emission sometimes appears significantly 
displaced from the location of the ultraviolet sources, at a resolution of 38" (FWHM) or 
about 400 pc. In the outer disk, dust may be removed away from star-forming sites more 
efficiently than in the inner disk, either because the volume density of the interstellar medium 
is decreased or because the exciting clusters are older on average and have injected more 
non-thermal energy into the interstellar medium. 

There are several examples of very porous star-forming sites at distances of 3-4' (~ 
2 kpc) from the center, from which both the hydrogen gas and the dust have been removed 
(such as D118 and D137) ; as they belong to the optically-bright part of the disk, where 



-25 - 



Table 4. 24 /im sources with weak optical and no ultraviolet counterpart. 



RA 


DEC (J2000) 


identification 


(h m s) 


(° ' ") 




00 54 41.7 


-37 40 22 


ultracompact H II region ? 


00 54 50.4 


-37 38 51 


background galaxy ? 


00 54 50.7 


-37 41 55 


ultracompact H II region ? 


00 54 53.3 


-37 43 13 


background galaxy ? 


00 54 59.8 


-37 39 27 


ultracompact H II region ? 


00 55 01.7 


-37 39 32 


ultracompact H II region ? 


00 55 05.9 


-37 42 14 


D125 - ultracompact H II region ? 
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presumably the molecular gas density remains high, i.e. above the threshold for gravitational 
instabilities discussed by Kennicutt (1989), these individual H II complexes are unable to 
affect significantly the global UV/IR radial profile. However, such regions may become 
dominant in the outer disk. 

Most of the Ha emission in the outer disk is in the form of diffuse, very large shells, 
which supports either one of the factors invoked above for dust removal: either the gas is very 
rapidly pushed away, during the lifetime of ionizing stars, because the interstellar medium 
is more tenuous, or the exciting stellar populations are older on average, so that they had 
more time to blow the gas away by their winds and supernova explosions. Good supporting 
examples of the hypothesis of aged stellar populations are the clusters D115 and D147, which 
are bright in the FUV, but seem to contain slightly evolved stellar populations (of the order 
of 5Myr, as opposed to 3Myr for most clusters). Since it is not possible to measure and fit 
the SEDs of a reasonably large number of clusters associated with giant shells in the outer 
disk, however, we cannot constrain the importance of this scenario relatively to the others 
(but see also Figure 2). 

From population synthesis fitting of individual young stellar clusters or small groups of 
clusters, we find that the stellar attenuation law varies significantly on the explored spatial 
scales of ~ 100 pc, and we argue that this empirical result likely corresponds to real variations 
in the extinction law as well. Just as the UV/IR ratio does not depend only on age and 
metallicity, there is ample evidence that the apparent correlation between the shape of the 
extinction law and the metallicity or the mode of star formation (active versus quiescent) 
is only a secondary effect of a more direct cause. We argue that the strength of the 2175 A 
bump and the steepness of the far-ultraviolet rise, whatever the responsible dust species, are 
directly related to the degree of compactness of the star-forming regions, and to the column 
density of interstellar material in front of the clusters. While most of the young stellar 
clusters considered here have similar ages and have not produced supernova explosions yet, 
the only obvious difference between them is that some appear to have formed in regions 
deeply affected by previous generations of stars. These are associated with diffuse, shell-like 
H II regions of large dimensions. If we divide star-forming complexes in two categories, 
based on the degree of compactness of the Ha emission, we find that the most compact 
regions follow either the Milky Way or the average LMC extinction law, whereas the most 
diffuse regions follow an extinction law with a decreased 2175 A bump and/or a steeper far- 
ultraviolet rise. While we cannot ascribe these features to a particular dust species, we show, 
using the near-infrared dust component extracted from our SED fitting procedure, thanks to 
the favorable angular resolution, that the distinction between the two types of morphologies 
and the two types of extinction laws is related to the column density of near-infrared emitting 
dust lying in the line of sight of the clusters. 
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We can hypothesize that the main difference in dust nature causing the change of ex- 
tinction law between the compact and diffuse regions is related to a changing proportion 
of dust grains very close to the H II region, which have been processed by hard ultraviolet 
radiation from the young stars, with respect to grains lying further away and which have not 
been processed in the same way. This photo-processing is likely occurring very fast, so that 
the observed dust properties of compact H II regions should be systematically different from 
the dust properties of diffuse regions where dust lies further away from the exciting stars. 
According to Boulanger et al. (1990), who discussed the ubiquitous and large variations in 
mid-infrared to far-infrared colors within and among nearby molecular clouds, such varia- 
tions cannot be explained only by changes in the radiation field, but require actual variations 
in the abundance of small grains with respect to big grains. They proposed that these abun- 
dance variations are caused by photo-processing, on time scales comparable to those of the 
mixing of gas by turbulent motions in clouds. Habart et al. (2001) and Abergel et al. (2002) 
also discussed abundance variations of small grains and aromatic hydrocarbons in resolved 
photodissociation regions and molecular clouds. This speculation would be consistent with 
the fact that the dust emitting in the 3-5 /im range is likely excited inside or very close to 
H II regions (Helou et al. 2004), and that we have seen here that its excitation has a close 
correspondence to the shape of the extinction law. 

A. Notes on individual regions (including the clusters for which the fits are 

shown in Figure 7) 

D53a: Compact, associated with very hot dust at 24 /im and 3-5 /im. Since the emis- 
sion is confused by nearby clusters in the GALEX bands, we had to subtract a gaussian 
representation of these contaminating clusters, and apply a large aperture correction to the 
FUV and NUV fluxes, because of the small size of the aperture. 

D53b: Compact, associated with hot dust at 24 /im and 3-5 /im. Two distinct clusters 
can be distinguished, but are well fitted with a single stellar population. 

D53c: Two clusters, one with a very red SED. The morphology in the various bands 
clearly indicates that interstellar emission is negligible at 3-5 /im. The fits included a second 
population but were not constrained to reproduce the 3-5 /im fluxes. The modelling confirms, 
however, the stellar nature of the 3-5 /im emission. The ionized gas is diffuse and forms an 
extended fragmentary shell of size ~ 200 pc. 
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D118: Complex region with several clusters, which can be fitted by a single population. 
The ionized gas and the dust are distributed in both compact regions and a giant shell. 

D119: Same remarks as for D118. 

D137b: This complex region stands out as very peculiar: the ionized gas and dust seem 
to form two walls on each side of the brightest cluster, separated by ~ 80 pc in projected 
distance. Star formation has thus had dramatic effects on the shaping of the surrounding 
interstellar medium. The Ha and dust peaks are coincident with weak UV sources, presum- 
ably because of extinction effects, while the cluster dominating the UV emission is located in 
a zone nearly void of interstellar emission. Since the complex contains redder stellar sources 
besides these clusters, we allow for a second population in the fits. This increases dramati- 
cally the degeneracy of the younger population parameters. 

D84: Compact region, associated with hot dust at 24 /im and 3-5 /im. 
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Fig. I. — (a) Radial profiles of the FUV and infrared surface brightnesses, measured in 
elliptical annulli oriented along the galaxy major axis. The black line corresponds to the 
FUV band and the grey line connecting the diamonds to the total infrared emission estimated 
from the 24, 70 and 160 /zm maps. The curves are normalized by their maximum value, (b) 
Radial profile of the FUV to infrared energy ratio, before (star symbols) and after (diamond 
symbols) correction of the FUV map for extinction. The energy emitted in the FUV band 
was computed with a filter equivalent width of 268 A. The linear dashed lines represent the 
opposite of the metallicity gradient derived by Deharveng et al. (1988) (with its uncertainty 
range). The big grey diamonds are individual FUV sources, measured at the same angular 
resolution of 38", within 76" apertures ; their size is proportional to the square root of the 
FUV flux corrected for extinction, using the results from Figure 4. 
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Fig. 2. — Radial profile of the FUV to NUV energy ratio, computed in the same way and 
at the same angular resolution as the FUV/TIR profile. Since the differential extinction 
between the two bands is very small, the extinction-corrected profile, if plotted, would be 
almost indistinguishable. We used filter equivalent widths of 268 A and 732 A, respectively. 
The flux density ratio is higher by a factor 1.22, or 0.087 dex. The error bars represent the 
3<t uncertainties on the ratio. 




Fig. 3. — Ha mosaic, with circles around the H II complexes discussed in Sect. 4 superposed. 
Their diameter is here arbitrarily fixed at 1', for clarity. The horizontal axis represents 
the right ascension and the vertical axis the declination, in the J2000 equatorial coordinate 
system. 
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Fig. 4. — Extinction in the FUV band (1516 A) derived from an energy balance argument: 
the power absorbed between 0.1 and 3/zm is set equal to the total infrared power. The 
assumed extinction law is that of Cardelli et al. (1989), but results are very similar with the 
other laws considered in this paper. The superposed circles represent the locations of the 
H II complexes discussed in Sect. 4. The coordinate axes are as in Figure 3. 




Fig. 5. — H I flux map from Puche et al. (1990), with the same coordinate axes and the 
same superposed circles as in Figures 3 and 4. 
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Fig. 6. — Maps of the complexes in U (preferred to FUV and NUV because of the more advantageous 
angular resolution), in Ha, and at 8^im. The horizontal axis represents the right ascension and the vertical 
axis the declination, in the J2000 equatorial coordinate system. The scale varies from one complex to another. 
Ha contours are superposed onto the U-band and 8 /im images. The intensity scale is linear for the 8 ^tm 
images, square root for U and Ha. Figure abridged for astro-ph submission. 





Fig. 7. — Fits of a few stellar cluster SEDs. The line connecting diamond symbols with 
error bars represents the observed SED after removing the Ha and [O II] contributions from 
the R and U bands ; the connected symbols without error bars represent the total (line + 
stellar) R and U-band fluxes. The dashed line is the best fit with the relevant extinction 
law, and the grey shading shows the la range of solutions. When present, two additional 
dashed lines indicate that two populations of different ages were fitted and represent their 
separate SEDs. 
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g. 7. — (continued) 
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Fig. 8. — Confidence interval plots for the fits shown in Figure 7. The areas of increasing grey intensity 
represent the 3a, 2a and la intervals without any constraint on the ionizing photon flux error (see text), 
and the black rectangle indicates the best fit. Due to the discrete nature of the extinction and age grids, 
individual solutions with error intervals (±0.05 mag, ±0.5 Myr and ±5% for the mass) are plotted rather 
than contours. The hatched areas correspond to solutions predicting an ionizing photon flux within 30% 
of the value estimated from the Ha and H/3 maps. No constraint on the FUV/NUV color is implied in x 2 
values. 
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Fig. 9. — Left: Relationship between the nebular extinction at 6563 A, estimated from the 
Ha/H/3 decrement, and the stellar extinction at 1516 A derived from the SED fits. The 
dashed lines represent the prediction of the ratio of stellar to nebular extinction A(R) = 
0.5A(Ha) from Calzetti et al. (1994), extrapolated to the FUV band using the four extinction 
laws considered in this paper. Right: Ratio of the nebular extinction at 6563 A to the 
stellar extinction in the R band, removing the ambiguity on the extinction law which was 
appropriately chosen for each cluster, according to the fit results. The dashed line shows 
again the prediction from Calzetti et al. (1994). The two points at A(Ha) = 0.28 and 
A(H«)/A(R) ~ 4.4 are D6 and D7-8, which lie outside the H/3 map and for which we 
adopted the Ha/H/3 decrement of Webster & Smith (1983) ; the fitted stellar extinction of 
D7-8, being zero, was replaced by its upper limit (Table 3). 
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Fig. 10. — Extinction and attenuation laws predicted by the DIRTY model, with an opacity 
Ty = 0.25 (corresponding to A 15ie = 0.7 using the Cardelli extinction law) for the Milky 
Way (thick black curves), the LMC2 region (thick grey curves) and the SMC (thin black 
curves). The extinction laws are shown as continuous lines ; the attenuation laws in the 
homogeneous dust distribution case as dashed lines ; and the attenuation laws in the clumpy 
dust distribution case as three- dots- dash lines. The approximate full widths at half maximum 
of the Galex transmission curves are indicated. 
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Fig. 11. — Predictions of the DIRTY model for the observed FUV/NUV flux density ratio of a 
3 Myr-old cluster, as a function of the extinction at 1516 A. The continuous lines represent the 
pure extinction laws folded in the DIRTY model: Milky Way, LMC2 (called 30 Doradus here) 
and SMC. The dashed lines represent the corresponding attenuation laws in the homogeneous 
case, and the three-dots-dash lines the attenuation laws in the clumpy case. The vertical 
double bars show the total range of FUV/NUV colors derived from each extinction law at 
^4i5i6 — 2.2, for easier identification. On each panel, the observed FUV/NUV ratios of the 
clusters well fitted with each considered extinction law (Cardelli, LMC and 30 Doradus) are 
overplotted. The clusters with unconstrained extinction laws (well fitted with all three laws) 
are omitted ; some clusters are in common for the Cardelli and LMC laws (see Table 3). The 
error bars account for the color variations expected when the age of the clusters is allowed 
to vary from 1 to 5 Myr. The maximal range of color variations, including these error bars, 
is indicated by a grey rectangle. 
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Fig. 12. — Relation between the Ha/R flux ratio, measured inside the stellar "aperture", 
and the inverse projected surface area of the entire H II region. The area was summed inside 
an aperture chosen to encompass as much of the H II region as possible without including 
nearby H II regions, and using only the pixels above 6a. 
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Fig. 13. — Compactness of the regions (as quantified by the Ha/R flux ratio) and extinction 
laws that are able to reproduce the observed spectral energy distributions from the FUV to 
the I band (see text). The label "not constr." means that the extinction law shape is not 
constrained by the stellar SED fits. 
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Fig. 14. — Ratio of the estimated interstellar component flux at 3.6-4.5 //m to the extinction- 
corrected FUV flux (chosen as normalization because it should be closely correlated with 
the heating intensity of the considered dust species), and applicable extinction laws. 
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Fig. 15. — Correlation between the interstellar component at 3.6 /xm and the aromatic band 
emission at 8 jum, measured within the stellar apertures, at the 8 /xm angular resolution. The 
dashed lines show brackets on the interstellar F 3 . 6 /F 8 flux ratio derived by Lu (2004). The 
relation with the 24 /xm emission is apparently less tight, but we are severely limited by the 
degraded angular resolution at this wavelength, about twice as large as at 8 /xm. 
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Fig. 16. — Correlation between the emission at 8 and 24 /im and the emission at 70 and 
160 /im, at the angular resolution of the 160 /im image. The solid lines indicate the least 
absolute deviation fits, and the dashed lines have a slope of 1. 



